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Abstract 

The structure and magnitude of the electromagnetic field produced by a rotating 
. accretion disk around a black hole were determined. The disk matter is assumed to 

O ' be a magnetized plasma with a frozenin poloidal magnetic field. The vacuum 

^ ■ approximation is used outside the disk. 

S3 ■ 
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^ INTRODUCTION 

in ; 

^^lecently, various models of particle acceleration near supermassive black holes (SMBHs) in galactic 
i — nuclei and near stellar-mass black holes (BHs) in the Galaxy have been widely discussed in connection 
^vith the studies of synchrotron radiation and inverse Compton scattering from narrow-beam jets 
Q)bserved over a wide spectral range, from radio to gamma rays. Nevertheless, the particularly high 
^ingular resolution provided by radio interferometers does not allow the central part of a quasar to 
Cbe distinguished, suggesting that the jet width is extremely small (comparable to the gravitational 
Qadius). Previously (Shatskiy 2001), the mechanism of Blandford and Znajek (1977) and Blandford 
^3{2001) for electric-. eld generation through the interaction of the magnetic .eld from a ring current 
c^vith the gravimagnetic .eld (GMF) of a Kerr BH located on the common axis with the ring current 
;>was used as the model of particle acceleration. Bisnovatyi-Kogan and Blinnikov (1972) and Shatskiy 
►v*hd Kardashev (2002) considered the mechanism of Deutsch (1955) or Goldreich and Julian (1969) 
$_for electric-field generation by the unipolar induction produced by the axial rotation of an accretion 
with a frozen-in magnetic field. Here, we suggest a mechanism that combines a unipolar inductor 
and strong gravitational SMBHs effects. Naturally, this mechanism is closer to the actual processes 
that take place in quasars. In contrast to the mechanism from Blandford and Znajek (1977), Beskin 
et al. (1992), and Beskin (1997), the mechanisms of Bisnovatyi-Kogan and Blinnikov (1972) used 
previously (Shatskiy 2001; Shatskiy and Kardashev 2002) operate in the vacuum approximation; the 
validity criterion for the latter is the condition of Goldreich and Julian (1969) for the number density 
of free charges: 

n e < |(nH)|/(27rce). (1) 

Here, Q is the plasma angular velocity, H is the characteristic magnetic field, c is the speed of light, 
and e is the elementary charge. Shatskiy and Kardashev (2002) showed that condition (1) could be 
satisfied near a BH, because there are no stable orbits for the particles closer than three gravitational 
radii (in a Schwarzschild field). The matter for which the vacuum approximation breaks down must 
be in the accretion disk, which, because of the effect of Bardeen and Petterson (1975), must be 
located in the equatorial plane of a rotating BH. We do not consider models in which the vacuum 
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approximation breaks down in the entire space outside a BH (models with the magnetohydrodynamic 
approximation). These were considered in detail, for example, in the review articles by Beskin et al. 
(1992) and Beskin (1997). Here, we determine the energy of the charged particles accelerated by the 
unipolar mechanism as well as the configurations of the electromagnetic field and the acceleration 
region. 



2 CONSTRUCTING THE MODEL 

Consider a Schwarzschild BH surrounded by an equatorial accretion disk. Let the disk have the 
characteristic size R and width 2a (see the figure). Because of the Bardeen-Petterson effect, the 
disk thickness can be disregarded. We use the following notation: M is the mass of the central 
body, R g = 2M is the Schwarzschild radius, 2 m is the mass of the test particle, Uj is its 4-velocity, 
Fij = diAj — djAi is the electromagnetic-field (EMF) tensor, Aj is the EMF potential, and Y l - k 
are the Christoffel symbols. Let us write the Schwarzschild metric and its determinant in spherical 
coordinate: 

ds 2 = (1 - r g /r)dt 2 - (1 - r g /r)- l dr 2 - r 2 d9 2 - r 2 sin 2 6d<p 2 , 

(2) 

g = — r 4 sin 2 9 . 

In general relativity, the following quantity for a charged particle that moves in stationary fields is 
conserved: 

e = m(u - 1) + eA , (3) 

which matches the particle energy in the nonrelativistic case. To prove this, it will suffice to consider 
the equation of motion for a charged particle in general relativity (see Landau and Lifshitz 1988): 

m lls = mukuiVlik + eykF ^ ( 4 ) 

where ds is the element of the particle proper time [see formulas (2)]. After transformation, this 
expression reduces to 

— {mui + eAi) = ^u l u k dig k i + eu k diA k . (5) 
ds I 

The conservation of energy . throughout the particle evolution follows for % = in stationary fields. 

Since the magnetic field is frozen into the disk, its distribution inside the disk is determined only 
by the initial conditions of the problem. These conditions depend on the accretion- disk formation 
mechanism. If the disk is assumed to have been formed through the destruction of a star by BH 
tidal forces, then the magnetic field of this star in the disk will preserve its direction. This field can 
have the profile shown in the figure. 

In the frame of reference comoving with the accretion disk, there is no electric field inside the disk 
because of its conductivity. Therefore, in a fixed frame of reference (with respect to distant stars), 
an electric field is induced by disk rotation inside the disk. Let the plasma in the disk rotate at an 
angular velocity Q relative to distant stars. The transformation of coordinates to a rotating frame is 
then: 3 

d^ = dx' k [5l + ^l5l5t}. (6) 



2 Below, we use the system of units in which the speed of light and the gravitational constant are equal to unity: c 
= 1, G = 1. 

3 Unlcss otherwise specified, Xi — t, r, 6, <p; x a = R,9 (Greek and Roman indices). 
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Because of axial symmetry, only the following EMF potential components are nonzero: Aq, the 
electric- field potential, and A v , the magnetic-field potential. According to (6) (see Landau and 
Lifshitz 1988), the EMF components transform as 

A' = A + VtA^, A'^ = A v , F' aQ = F a0 + ftF aip , F' aif = F atp . (7) 

Since F a0 = in plasma, we have inside the disk 

F a0 = -QF a(p , A = const - QA^. (8) 

On the disk surface, continuous boundary conditions exist for the tangential electric-field components 
and for the normal magnetic-field components. Outside the disk, there are no field sources by the 
definition of the model. Thus, determining the EMF reduces to solving the Laplace equation in the 
spacetime curved by gravity with the specified boundary conditions on the disk surface and on the 
BH horizon. The boundary conditions for the EMF tensor on the BH horizon were found previously 
(Shatskiy 2001): 

lim F oe oc g o -> , lim F TLp oc goo -> . (9) 

r— >r g r—>r g 

In turn, the boundary conditions on the disk for the magnetic and electric fields are determined 
solely by the magnetic- field distribution inside it. The specific form of this distribution is not that 
important for the solution of the problem. This is because at distances from the disk to the point of 
observation much larger than the disk thickness, the dipole field (the monopole field must be absent, 
because the total disk charge is zero) mainly contributes to the electric field of the disk element 
within position angles between tp and tp + dtp when the field is expanded in multipoles. In this case, 
the total electric field obtained by integrating over the angle tp has a quadrupole nature: 

lim A = const ■ (1 — 3 cos 2 9)/r 3 . 

(r/r g )^oo 



3 THE MAXWELL EQUATIONS 

The Maxwell equations for the EMF in general relativity are 

-±=d l {^F' lk )=An 3 k . (10) 

In the disk, the 4-vector of the current j k can be determined from a given magnetic field. Let us 
introduce the physical components of the EMF vectors, their analogs in Euclidean space: 4 

E a = -^ 0V VVa H a = -e Q ^F wV /|^V1, J"=/\fe. (11) 

Here, e al3y = e a/3l is the Levi-Civita symbol. This form of the EMF physical components was chosen 
in order that Eq. (10) correspond to the classical Maxwell equations in Euclidean space: 

div E = 47rj°, rot H = 4ttJ. (12) 

4 We denote them by a hat. 
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4 THE ELECTROMAGNETIC FIELD NEAR A SMBH 

The magnetic field of an accretion disk around a Schwarzschild BH was determined by Tomimatsu 
and Takahashi (2001). The electric field of the disk element within position angles between ip and 
ip + dip can be represented as the field from two charges: +<?ff and — q& , located inside the disk, 
at the system equator, and at distances +a and —a from its center (r = R, 9. = 7r/2), respectively 
5 . As a result, we obtain an electric dipole 2qadf/7i in the disk element between ip and ip + dip. In 
Euclidean space, the radial electric field at the disk center can be obtained by integrating over the 
angle ip; for a « R, it is E = E^ r=Re=n ^ 2 ^ = —2q/(nRa). The corresponding magnetic field (which 
is responsible for the emergence of the electric field) can be found from the electric field. Note that 
nothing forbids the frozenin magnetic field in the disk to have precisely such a profile (see the figure). 
In the figure, the accretion disk is located at r ~ QM. The contradictions related to the existence of 
stable orbits in this region can be removed by the following reasoning: 

1 For the Kerr metrics, the nearest stable orbit is at radius r = M (see, e.g., Landau and Lifshitz 
1988). 

2 Even if the orbit is not stable, it is spiral and goes under the horizon, while a new orbit can 
arrive in place of it. Thus, the pattern is quasi-stationary. According to (11) and (8), the quantity 
q can be expressed in terms of the magnetic field at the disk center 6 as follows: 



7T 



—TrRaEo/2 = -ttR 2 aH . 



(13) 



In physical coordinates, the field of a point charge e near a BH was presented by Thorne et al. (1998). 
It was obtained in a closed form by Linet (1976): 



A = 



e 

Rr 



M | {R-M)(r-M)-M'h 



E r 



E e 



1 - 



R—M+Mt 



D 



+ 



r[(r-M)(R-M)-M 2 t][r-M-(R-M)t] ~) 
D3 J 



e(fl-2M) v /l-2M/r a 
O t , 



(14) 



D 2 = (r - Mf + (R- Mf — M 2 — 2(r - M)(R - M)t + MH 2 . 

where t is the cosine of the angle between the directions of the point charge and the point of 
observation of the field from the BH center. The electric field of a charged ring at the BH equator 
was found by Bicak and Dvorak (1996) in the form of a series. Here, this field is found in a quadrature 
form. To this end, we make the following substitutions: t — > sin 9 cosy?, e — > Q^- and integrate over 
the angle from — n to +ir. The model electric-field potential A^ 1 is a superposition of the fields 
from two charged rings at radii R + a and R — a and with charges +q and —q, respectively. The 
quadrature obtained can be expressed in terms of incomplete elliptic integrals. Since this quadrature 
is cumbersome, it makes no sense to write it here. Instead, we give an expression more useful for 
practical calculations, an expansion of this quadrature in terms of a. We retain only the first term 



5 Naturally, there are no free charges in the disk; these were introduced for the convenience of representing the 
dipole field outside the disk. 

6 H = -He( r =R,8=ir/2)l \J 1 - r g/R- 
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of the series (because a « R is small): 7 

+7T +7T 



Al ot = j 2ad R A dip= J Aa(d r A ) (R ^ r) dip = J AaE[ R ^ r) d V . 



Substituting expression (14) here finally yields 



Atot _ 2aq } A ( M(D—r+M—Mt) , R[(r-M)(R-M)-M 2 t}[R-M-(r-M)t] \ 

^0 - nrR 2 l ai P\ D + D 3 I • {i-O) 





On the Q axis, the integration over ip is simple and it is easy to see that the only maximum of the 
potential Aff* (outside the BH) is on the horizon (r = 2M): 

K«(r = 2tf,« = 0)L„ = 22^. (16) 

Consider the distribution of energies e(r) in an ensemble of test particles at rest on the Q axis 
[see (3)]. This distribution has a weak maximum near the horizon and then slowly falls o. to zero at 
distances r » r g . Subsequently, the energy variation can be virtually disregarded. 



5 DISCUSSION 

The necessary condition for particle escape to infinity from rest is positiveness of the force rn^- = 
that acts on the test particle at the starting point. This requires that the charge of the test 
particle e have the sign opposite to that of q: 

- sign (eg) = sign (ef2E (e=0) ) = - sign (ef2H ) = 1. (17) 

The force that acts on an particle is proportional to the particle energy gradient (see Eq. (4)]: 

du r i 

m— = -{d r e)- y / l-2M/r. (18) 

(JjO 

In the main order in a = mM/(eq) « 1, the energy maximum is almost equal to the electric 
component of the particle energy on the horizon: 8 

2Ma \ , 9 , „ , 

-1 -mc 2 ~ 0.1 -mc 2 /a (19) 



,a B? — a 2 

and the point at which this maximum is reached is at the following distance from the horizon: 

2 (R 2 - a 2 ) 2 [(R - M) 2 - a 2 } 4 3 2 

r max r g -r 9 -a ■ 16M6(l2 ^ R2 + M 2 _ 4MR + a 2y ~ r g ■ 10 ■ a . (20) 

In conclusion, several more words can be said about the same model in a Kerr field. BH rotation gives 
rise to a GMF that interacts with the EMF of the disk and changes its components in magnitude and 
direction. In the linear approximation in BH angular velocity, the GMF gives an additive component 

7 Here, we make use of the symmetry of the potential A in variables R and r and use their change (the subscript 

It 55 \ 

(fi~r) )■ 

8 At a~ r g ~ R/3 [see (16)]. 
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Figure 1: Accretion disk around black hole, which center at radius R = 3, and which width a = 1, 
with the frozen-in magnetic field — firm line. Electrical field (generated by magnetic) — dotted line. 
Z and p are expressed in fractions of r g = 2M. 



to expression (19) for the maximum energy of the charged particle accelerated by an electric field. 
This component was determined previously (Shatskiy 2001; it is convenient to represent it here as 

Here, Q g is the angular velocity of the BH horizon (the falling test particles are drawn into rotation 
by the BH GMF).We see from (19) that at 

the particle energy 9 accelerated by SMBHs can reach values larger than 10 20 eV. 
9 Here, m e is the electron mass. 
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6 CONCLUSIONS 



1 The model described above yields even a higher energy of the accelerated particles than does the 
Blandford-Znajek model. 

2 The energy excess can be expended on the radiation reaction and on collisions with particles of 
the rarefied plasma near the axis. 

3 The derived EMF configuration can be used to numerically calculate the dynamics of the 
charged particles far from the axis. It will make it possible to compare theoretical conclusions with 
observational data. (4) Far from the axis, a force-free field with the approximation of magnetohy- 
drodynamic models can give a large contribution to the EMF amplitude. Therefore, in numerical 
calculations, the contributions to the EMF configuration from different models should be taken into 
account with different weights. 



This work was supported by the Russian Foundation for Basic Research (project nos. 01-02-16812, 
00-15-96698, 01-02-17829). 

7 ACKNOWLEDGMENTS 

It wish to thank N.S. Kardashev, R.F. Polishchuk, V.N. Lukash, B.V. Romberg, Yu.Yu. Rovalev, 
the remaining sta. of the theoretical departments at the Astrospace Center and the Lebedev Physical 
Institute, and participants of workshops for an active participation in preparing the paper and for 
important remarks. 

8 REFERENCES 

[I] . J.M.Bardeen and J. A. Petterson,Astrophys. J. Lett. 195, L65, 1975. 
[2]. V. S. Beskin, Usp. Fiz. Nauk 167, 689, 1997. 

[3]. V. S. Beskin, Ya. N. Istomin, and V. I. Par'ev, Astron. Zh. 69, 1258, 1992. 
[4]. G. S. Bisnovatyi-Rogan and S. I. Blinnikov, Astrophys. Space Sci. 19, 119, 1972. 
[5]. J. Bicak and L. Dvorak, General Relativity and Gravitation 7, 959, 1996. 
[6]. R. Blandford and R. Znajek, MNRAS 179, 433, 1977. 

[7]. R. Blandford, Galaxies and Their Constituents at the Highest Angular Resolution. Proc. 
IAI Symp. 205. Ed. R. T. Schilizzi (San Francisco, ASP, 2001), p. 10. 
[8]. P. Goldreich andW. H. Julian, Astrophys. J. 157, 869, 1969. 
[9]. J. Deutsch, Ann. D'Astrophys. 1, 1, 1955. 

[10]. L. D. Landau and E.M. Lifshitz, Field Theory (Nauka, Moscow, 1988). 

[II] . B. Linet, J. Phys. A 9, 1081, 1976. 

[12]. A. Tomimatsu and M. Takahashi, Astrophys. J. 552, 710, 2001. 

[13]. R. Thorne, P. Price, and D. MacDonald (Eds.), Black Holes. A Membrane Approach 
(Mir,Moscow, 1998). 

[14]. A. A. Shatskiy, Zh. Eksp. Teor. Fiz. 93, 920, 2001, (gr-qc/0202068). 

[15]. A. A. Shatskiy and N. S. Rardashev, Astron. Zh. 46, 639, 2002, (astro-ph/0209465). 



Translated by G. Rudnitskii 



7 



